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EXECUTIVE  SUMMARY 


The  funnel-and-gate  system  is  an  in  situ  technique  for  the  restoration  of  contaminated 
aquifers.  It  consists  of  low  to  impermeable  barriers  used  to  divert  the  groundwater  flow  through 
a  controlled  area.  To  gain  insight  to  the  flow  field  dynamics  of  such  systems,  a  series  of  three- 
dimensional  calculations  was  performed.  In  the  calculations  the  width  of  the  gate  and  the  depth 
of  the  barriers  were  varied.  The  individual  calculations  were  designed  using  the  Plackett- 
Burman  statistical  experimental  design  method.  The  results  of  the  calculations  were  used  to 
develop  a  statistical  model  of  the  fuimel-and-gate  system  (reported  separately).  The 
computations  show  little  effect  from  varying  the  angle  of  the  walls  to  the  flow  field  in  capturing 
the  flow.  The  capture  zone  is  also  imaffected  by  increasing  the  hydraulic  head.  The  higher  the 
conductivity  at  the  gate,  the  greater  the  flow  and  also  the  capture  zone. 
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SECTION  I 


INTRODUCTION 


OBJECTIVE 

The  purpose  of  this  work  was  to  examine,  in  3D,  groundwater  flow  in  the  vicinity  of  a 
partially  penetrating  funnel-and-gate  system. 

BACKGROUND 

The  funnel-and-gate  system  is  an  in  situ  technique  for  the  restoration  of  contaminated 
aquifers.  It  consists  of  low-permeability  barriers  or  cutoff  walls  (funnels)  which  divert  the  flow 
of  groundwater  towards  and  through  the  gate.  The  gate  consists  of  one  or  a  sequence  of  reactive 
materials  which  can  either  degrade  or  retard  the  groundwater  contaminants.  A  funnel-and-gate 
system  can  be  either  fully  penetrating  or  partially  penetrating.  A  fully  penetrating  funnel-and- 
gate  system  implies  that  it  is  installed  from  the  ground  surface  and  anchored  at  some  depth  below 
ground  surface  into  a  low  permeability  material  (e.g.,  clay  layer).  A  partially  penetrating  system 
is  one  that  is  not  anchored  into  a  low  permeability  material,  and  thus  groimdwater  can  flow  freely 
vmder  the  system.  Figure  1  shows  schematic  diagrams  of  a  partially  penetrating  funnel-and-gate 
system  in  plan  view  (Figure  la)  and  cross  section  (Figure  lb). 

To  gain  insight  into  the  effectiveness  of  funnel-and-gate  systems,  a  numerical  modelling 
study  is  required.  Preliminary  work  on  modelling  funnel-and-gate  systems  includes  that  by  Starr 
and  Cherry  (1994)  who  performed  a  detailed  sensitivity  analysis  in  2D  plan  view,  and  Shikaze  et 
al.  (1995)  who  presented  results  in  3D.  Both  of  these  studies  simulated  steady-state  groundwater 
flow  in  the  vicinity  of  a  funnel-and-gate  system.  Christenson  and  Hatfield  (1994)  developed  an 
analytical  model  to  examine  funnel-and-gate  systems.  They  used  Dupuit-Forcheimer 
assumptions  and  developed  plots  showing  the  relationships  between  various  geometrical 
parameters  of  a  funnel-and-gate  system  in  2D  areal  view.  It  is  worth  noting  here  that  all  previous 
3D  work  assumes  that  the  funnel-and-gate  system  is  fully  penetrating.  The  key  difference  for  the 
2D  simulations  is  that  flow  under  the  cutoff  walls  can  be  accounted  for.  In  the  previously 
mentioned  2D  studies,  partially  penetrating  conditions  cannot  be  considered.  Shikaze  et  al. 
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funnel-and-gate  systems  in  3D.  However,  they  only  performed  a  limited  number  of  simulations 
and  only  examined  one  funnel-and-gate  geometry. 

Figure  2a  shows  a  perspective  view  of  a  typical  3D  domain.  In  this  case,  the  funnel-and- 
gate  system  is  partially  penetrating  and  the  dots  represent  the  beginning  of  a  particle  trace  in  a 
steady-state  groundwater  flow  field.  Figure  2b  shows  the  particles  again,  as  well  as  streamlines 
based  on  the  flow  field.  Here,  a  portion  of  the  original  particles  eventually  passes  through  the 
gate.  The  zone  of  particles  upgradient  that  ends  up  passing  through  the  gate  is  referred  to  in  this 
work  as  the  capture  zone,  and  is  shown  in  Figure  2b  as  a  black,  hatched  zone.  In  this  report, 
capture  zones  will  be  shown  in  a  manner  similar  to  that  shown  in  Figure  3.  This  view  is  looking 
towards  the  funnel-and-gate  system  in  the  direction  of  groundwater  flow.  From  this  viewpoint, 
the  capture  zone  can  easily  be  drawn  and  measured. 

SCOPE 

A  number  of  simulations  which  involved  varying  one  or  more  of  the  dimensionless 
parameters  that  are  used  to  define  the  geometry  of  a  funnel -and-gate  system  were  performed. 

The  simulation  results  presented  here  have  been  chosen  by  Applied  Research  Associates,  Inc. 
(ARA),  and  are  in  the  form  of  a  "run-matrix"  of  16  simulations.  In  addition,  three  simulations 
have  been  performed  to  be  used  by  ARA  to  compare  with  another  computer  program.  These 
three  simulations  are  described  in  detail  in  Section  III. 

Section  II  describes  the  numerical  model  that  was  used  for  this  work,  including  a 
discussion  of  the  governing  equations  and  the  incorporation  of  impermeable  cutoff  walls. 

Section  III  presents  the  simulations  that  are  to  be  used  by  ARA  for  comparison  to  another  model, 
and  Section  IV  contains  the  results  from  the  simulations  for  the  run-matrix,  as  well  as  a 
definition  of  boundary  conditions,  and  dimensionless  parameters. 
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SECTION  II 


I 


THEORETICAL  BACKGROUND 

The  model  that  was  used  for  this  project,  FRAC3DVS,  is  described  in  detail  by  Therrien 
(1992)  and  Therrien  and  Sudicky  (1995).  It  is  a  three-dimensional  numerical  model  which  can 
be  used  to  simulate  saturated  or  partially  saturated  groundwater  flow  and  solute  transport  in 
porous  or  fractured,  porous  geologic  media.  For  the  purposes  of  this  work,  we  used  the  model  to 
simulate  fully  fractured,  steady-state  groimdwater  flow  in  a  homogeneous  nonfractured  aquifer 
(except  for  Case  2  in  Section  III,  where  steady-state,  partially  saturated  conditions  were 
simulated). 

NUMERICAL  MODEL 


wemmg  bauations 


The  partial  differential  equation  that  governs  fully  saturated,  steady-state  groundwater 

flow  is: 

d^h  d^h  d^h 

+  K,—^  =  0  (1) 

dx  ^dy  dz 

where  K^,  Ky,  are  the  hydraulic  conductivities  in  the  x,  y,  and  z  directions,  respectively,  h  is 
the  hydraulic  head,  and  x,  y,  and  z  are  the  spatial  coordinates. 

The  specific  discharge  is  calculated  by  Darcy's  equation: 
q  =  -Ki  (2) 

where  i  is  the  hydraulic  gradient,  and  q  is  the  specific  discharge. 

One  simulation  in  Section  III  involves  partially  saturated  groundwater  flow.  For  partially 
saturated,  steady-state  groundwater  flow,  a  modified  form  of  Richards'  equation  is  used  to 
describe  groundwater  flow  (Therrien,  1992): 


d  (  5  (v|/  +  z)^ 


i,j  =  x,y,z 


where  Ky  is  the  saturated  hydraulic  conductivity  tensor,  =  Kw(^J  is  the  relative  permeability 
of  the  medium  as  a  function  of  water  saturation,  =  i|/(x„  t)  is  the  pressure  head,  and  z  is  the 
elevation  head.  The  saturation  of  water  is  related  to  water  content,  0,  by: 
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(4) 


where  Of  is  the  residual  water  content,  and  Og  is  the  porosity.  For  the  partially  saturated  example 
in  Section  III,  we  use  tabular  data  as  input  for  the  Sy^/(y/)  and  kry^(Syv)  relationships. 


Within  the  3D  computational  domain,  impermeable  cutoff  walls  are  implemented  as  2D 
planes.  This  is  done  by  adding  "false  nodes"  wherever  impermeable  nodes  are  desired.  As  a 
consequence,  wherever  an  impermeable  wall  exists,  there  are  two  nodes  existing  at  the  same 
spatial  location,  one  of  which  is  connected  to  an  element  on  one  side  of  the  wall,  while  the  other 
is  connected  to  the  element  on  the  opposite  side  of  the  wall.  The  result  is  to  essentially  break  the 
connection  between  two  adjacent  elements  thereby  mimicking  an  impermeable  wall  by 
emplacing  a  2D  impermeable  plane  within  the  3D  domain.  The  3D  block  elements  remain 
unchanged,  the  only  difference  being  new  elemental  incidences  for  elements  that  are  associated 
with  cutoff  walls. 
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SECTION  III 


THREE  COMPARISON  SIMULATIONS 

The  purpose  of  these  three  simulations  using  FRAC3DVS  is  to  present  contrasting  cases 
that  can  be  compared,  by  ARA,  to  results  from  another  program.  All  three  of  these  simulations 
represent  partially  penetrating  systems.  These  cases  are  (1)  fully  saturated,  steady-state 
groundwater  flow  with  cutoff  walls  perpendicular  to  the  direction  of  ambient  groundwater  flow, 
(2)  identical  to  case  (1)  except  partially  saturated  flow  conditions,  and  (3)  identical  to  case  (1) 
except  that  cutoff  walls  are  aligned  45  degrees  to  the  direction  of  ambient  groundwater  flow  such 
that  the  two  cutoff  walls  of  the  funnel-and-gate  system  form  a  "V"  shape.  The  computational 
domain  is  essentially  identical  to  that  used  in  the  simulations  from  the  Run  Matrix,  with  the  only 
difference  being  the  thickness  (i.e.,  in  the  z-direction)  of  the  domain.  The  three  simulations  are 
described  in  further  detail  below. 

CASE  1;  FULLY  SATURATED  FLOW;  CUTOFF  WALLS  PERPENDICULAR  TO  FLOW 

For  this  simulation,  the  parameters  and  boundary  conditions  are  listed  below: 

•  Domain:  100  x  50  x  10  meters  (in  x,  y,  z) 

•  Funnel-and-gate  depth  =  5  meters 

•  Gate  width  =  2  meters 

•  Cutoff  walls:  2-  x  10-meter  walls  (partially  penetrating) 

•  (1)  XY  plane  at  y  =  20  meters,  from  x  =  39  to  49  meters,  z  =  5  to  10  meters 

•  (2)  XZ  plane  at  y  =  20  meters,  from  x  =  5 1  to  61  meters,  z  =  5  to  10  meters 

•  Hydraulic  gradient  =  0.005  m/m  (flow  is  parallel  to  y-axis) 

•  ^aquifer  ~  10^  porosity  =  0.33 

•  ^gate  ~  ^aquifer 

•  Boundary  conditions: 

•  hydraulic  head  (hi)  =  10.25  meters  on  XZ  plane  at  y  =  50  meters,  from  x  =  0 
to  100  meters,  z  =  0  to  10  meters 

•  hydraulic  head  (h2)  =  1 0.00  meters  on  XZ  plane  at  y  =  0  meters,  from  x  =  0  to 
100  meters,  z  =  0  to  10  meters 
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•  all  other  domain  boundaries  are  impermeable 
•  number  of  nodes  is  approximately  1 80,000 

For  Case  1,  the  flow  rate  through  the  gate  was  determined  to  be  approximately  0.86 
mVday  (NOTE:  these  flow  calculations  were  determined  by  summing  the  specific  discharge,  q, 
of  all  3D  block  elements  at  the  front  face  of  the  gate,  and  multiplying  this  by  the  total  flow¬ 
through  area  of  the  gate.  A). 

The  resulting  capture  zone  for  Case  1  is  shown  in  Figure  4  and  is  parabolic  in  shape,  with 
decreasing  width  as  depth  increases.  The  parabolic  shape  is  a  result  of  the  fact  that  the  funnel- 
and-gate  system  is  a  hanging  system  and  groundwater  flow  is  permitted  to  pass  under  the  funnel- 
and-gate.  Table  1  lists  the  width  of  the  capture  zone  at  various  depths  below  the  top  of  the 
computational  domain. 


TABLE  1.  CAPTURE  ZONE  WIDTH  AT  VARIOUS  DEPTHS  -  CASE  1. 


Depth  (m) 

Capture  zone  width  (m) 

0.0 

11.5 

1.0 

9.5 

2.0 

7.3 

3.0 

5.1 

4.0 

2.6 

Other  results  from  Case  1  are  presented  in  Appendix  A,  including  2D  slices  of  velocity 
vectors  and  hydraulic  head  contours. 

CASE  2:  PARTIALLY  SATURATED  FLOW;  CUTOFF  WALLS  PERPENDICULAR  TO 
FLOW 

For  Case  2,  the  furmel-and-gate  geometry  has  been  set  up  to  be  identical  to  that  in  Case  1 
as  described  in  the  previous  section  (i.e.,  partially  penetrating  funnel-and-gate,  walls 
perpendicular  to  flow).  However,  Case  2  involves  partially  saturated  groundwater  flow,  and  as  a 
result  the  boundary  conditions  are  defined  differently,  as  presented  below: 

•  Boundary  conditions 

•  hi  =  8.25  meters 

•  h2  =  8.0  meters 

•  top  boundary  has  recharge  of  0.00055  m/day 
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Because  this  is  a  partially  saturated  simulation,  we  must  define  the  relationship  between 
pressure  head  and  saturation,  and  relative  permeability  and  saturation.  These  relationships  are 
tabulated  in  Table  2. 


TABLE  2.  CONSTITUTIVE  RELATIONSHIPS  FOR  UNSATURATED  ZONE 
PARAMETERS. 


Pressure  head  (m) 

Saturation 

Relative  k 

Saturation 

-10.0 

0.12 

0.0001 

0.20 

-3.0 

0.20 

0.01 

0.40 

-1.0 

0.35 

0.1 

0.63 

0.0 

1.0 

1.0 

1.0 

For  Case  2,  the  flow  rate  through  the  gate  is  approximately  0.78  mVday.  Figure  5  shows 
the  capture  zone  for  this  case  at  various  depths,  and  these  capture  zone  results  are  listed  in  Table 
3. 


TABLE  3.  CAPTURE  ZONE  WIDTH  AT  VARIOUS  DEPTHS  -  CASE  2. 


Depth  (m) 

Capture  zone  width  (m) 

0.0 

8.0 

1.0 

7.8 

2.0 

7.0 

3.0 

5.5 

4.0 

2.6 

Appendix  B  shows  other  figures  from  Case  2  (2D  plots  of  hydraulic  head  contours  and 
velocity  vectors). 

CASE  3;  FULLY  SATURATED  FLOW;  CUTOFF  WALLS  45  DEGREES  TO  FLOW 

Case  3  has  been  set  up  with  identical  boundary  conditions  to  Case  2  with  fully  saturated 
groundwater  flow.  However,  in  this  case,  the  two  cutoff  walls  are  oriented  45  degrees  to  the 
direction  of  groundwater  flow.  This  was  done  by  "stepping"  in  0.5-meter  increments  using 
orthogonal  cutoff  walls  (see  Figure  6).  As  before,  the  funnel-and-gate  system  penetrates  to  a 
depth  of  5  meters. 
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'  Figure  7  shows  the  resulting  capture  zone,  and  capture  zone  widths  are  listed  in  Table  4 
for  various  depths  below  surface. 


TABLE  4.  CAPTURE  ZONE  WIDTH  AT  VARIOUS  DEPTHS  -  CASE  3. 


Depth  (m) 

Capture  zone  width  (m) 

0.0 

13.5 

1.0 

11.1 

7.5 

2.2 

4.0 

0.0 

The  flow  rate  through  the  gate  is  approximately  0.75  mVday,  and  additional  results  for 
Case  3  are  shown  in  Appendix  C. 
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SECTION  IV 


SIMULATIONS  FOR  RUN  MATRIX 

COMPUTATIONAL  DOMAIN  AND  BOUNDARY  CONDITIONS 

The  computational  domain  for  simulations  for  the  run  matrix  has  been  set  up  to  be  100  by 
50  meters  in  the  x  and  y  directions,  respectively.  In  the  vertical  (z)  direction,  the  domain  was 
varied  between  10  and  18  meters  to  accommodate  the  desired  funnel-and-gate  geometry  for  each 
specific  simulation. 

For  all  simulations,  the  hydraulic  conductivity  and  porosity  of  the  aquifer  were  set  as  10"^ 
m/s  and  0.33,  respectively.  The  hydraulic  conductivity  of  the  gate  was  varied  from  a  low  of  10"^ 
m/s  to  a  high  of  200”^  m/s.  The  number  of  nodes  for  each  simulation  varied  such  that  an 
adequate  spatial  discretization  was  attained  in  the  vicinity  of  the  funnel-and-gate.  In  other  words, 
small  grid  blocks  were  used  near  the  furmel-and-gate,  and  larger  grid  blocks  were  used  further 
away  from  the  fimnel-and-gate.  The  range  in  the  number  of  nodes  was  between  about  164,000 
and  296,000. 

All  simulations  presented  in  this  section  involve  steady-state,  fully  saturated  groimdwater 
flow  in  a  3D  flow  domain  with  partially  penetrating  funnel-and-gate  systems.  Boundary 
conditions  have  been  set  up  such  that  the  ambient  flow  of  groundwater  is  parallel  to  the  y-axis 
(see  Figure  8).  This  is  done  by  setting  constant  head  boundary  conditions  along  the  end  XZ  faces 
(i.e.,  at  y  =  0  and  50  meters).  The  values  of  these  specified  head  faces  were  varied  in  order  to 
obtain  the  desired  hydraulic  gradient  for  each  simulation.  All  other  boundary  faces  were 
assumed  to  be  impermeable  boundaries.  • 

DIMENSIONLESS  PARAMETERS 

To  define  the  geometry  of  a  funnel-and-gate  system  as  well  as  relevant  aquifer 
parameters,  several  dimensionless  parameters  have  been  defined.  The  main  reason  for  using 
dimensionless  parameters  is  to  reduce  the  number  of  variables  in  the  system.  The  four 
dimensionless  parameters  are  as  follows:  (a)  the  ratio  of  the  hydraulic  conductivity  of  the  gate  to 
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that  of  the  aquifer, 
furinel-and-gate. 


K 


gate 


V  ^aquifer  J 


,  (2)  the  ratio  of  the  width  of  a  single  cutoff  wall  to  the  depth  of  the 


,  (3)  the  ratio  of  the  total  cutoff  wall  width  to  the  gate  width, 


2*  w 


(4)  the  hydraulic  gradient,  .  Figure  9  shows  the  definition  of  these  parameters  in  the 
context  of  a  typical  funnel-and-gate  system. 

RESULTS  FROM  RUN  MATRIX 


v^/y 


,  and 


Table  5  shows  the  16  simulations  as  shown  by  ARA.  These  simulations  were  chosen  on 
the  basis  of  minimum,  maximum,  and  mean  values  for  each  of  the  four  dimensionless 
parameters,  as  described  above.  Also  shown  in  Table  5  are  the  values  of  the  depth  of  the  funnel 
and  gate  (df),  as  well  as  the  width  of  one  cutoff  wall  (wf)  and  the  width  of  the  gate  (wg)  for  each 
simulation. 


TABLE  5.  SUMMARY  OF  PARAMETERS  USED  IN  THE  RUN  MATRIX. 


Run 

^gaie 

‘‘f 

2*Wf 

grad  h 

df{va) 

Wg  (m) 

wf{m) 

^aquifer 

1 

4 

2.0 

10 

5 

2 

10 

2 

1 

0.5 

5 

7 

1.4 

3.5 

3 

1 

0.5 

20 

0.010 

14 

0.7 

7 

4 

1 

4.0 

5 

0.010 

2.5 

4 

10 

5 

1 

4.0 

20 

5 

2 

20 

6 

0.5 

5 

0.010 

7 

1.4 

3.5 

7 

0.5 

20 

14 

7 

8 

4.0 

5 

2.5 

4 

10 

9 

4.0 

20 

5 

2 

20 

4 

4.0 

20 

0.005 

5 

2 

20 

11 

4.0 

10 

0.005 

2.5 

2 

10 

12 

1 

4.0 

5 

0.001 

2.5 

4 

10 

13 

4 

0.5 

10 

0.010 

10 

1 

5 

14 

4 

0.5 

20 

0.001 

14 

7 

15 

0.5 

10 

0.001 

10 

1 

5 

16 

1 

0.5 

5 

0.010 

7 

1.4 

3.5 
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Figures  10  to  25  show  the  capture  zones  for  these  16  simulations,  and  the  relative  capture  zones 
are  shown  in  Table  6. 

Qrel  (o’*  relative  flow  through  the  gate  is  defined  as  Qgate  divided  by  the  total  Q  through 
an  area  equal  to  the  fimnel-and-gate  area  (in  the  absence  of  the  funnel-and-gate).  Relative  depth 
is  defined  as  the  depth  of  measurement  of  the  capture  zone  width  divided  by  df,  and  the  relative 
capture  zone  width  is  defined  as  the  capture  zone  width  (in  meters)  divided  by  the  total  width  of 
the  funnel-and-gate  system  (which  is  defined  as  2  *  wf+  Wg). 

The  capture  zone  widths  at  relative  depths  of  0.0,  0.2,  0.4,  0.6  and  0.8  are  shown  on  the 
figures.  Because  the  measurement  of  the  width  of  a  capture  zone  is  approximate,  relative  capture 
zone  widths  are  accurate  only  to  within  about  0.02. 


TABLE  6.  RESULTS  FROM  THE  RUN  MATRIX. 


Run 

Qgate 

(—] 

Kday) 

Qrel 

Relative  capture  zone  width  at  relative  depth  = 

0.0 

0.2 

0.4 

0.6 

0.8 

1 

1.21 

0.45 

MSM 

WSBM 

2 

0.15 

■m 

HU 

0.49 

■H 

■H 

3 

3.03 

0.170 

0.43 

0.41 

0.41 

0.41 

4 

1.12 

0.216 

0.42 

0.35 

0.30 

0.23 

BB 

5 

0.19 

0.105 

0.29 

0.24 

0.18 

0.12 

■H 

6 

2.66 

0.524 

0.58 

0.58 

0.58 

0.58 

0.51 

7 

0.61 

0.343 

0.46 

0.46 

0.46 

0.46 

0.31 

8 

0.18 

0.347 

0.46 

0.40 

0.34 

0.28 

0.23 

9 

3.04 

0.168 

0.36 

0.27 

0.20 

0.14 

0.09 

10 

1.30 

0.144 

0.36 

0.27 

0.20 

0.14 

0.07 

11 

0.59 

0.249 

0.40 

0.35 

0.29 

0.23 

0.17 

12 

0.11 

0.217 

0.42 

0.35 

0.30 

0.23 

0.18 

13 

3.26 

0.343 

0.51 

0.51 

0.51 

0.51 

0.40 

14 

0.49 

0.275 

0.42 

0.42 

0.42 

0.42 

0.31 

15 

0.40 

0.418 

0.52 

0.52 

0.52 

0.52 

0.41 

16 

1.49 

0.298 

0.49 

0.49 

0.49 

0.49 

0.41 
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SECTION  V 


DISCUSSION 

The  effect  of  each  of  the  four  dimensionless  parameters  can  be  examined  by  comparing 
several  of  the  above  16  cases.  First,  Runs  2  and  16  can  be  compared  to  examine  the  effect  of 
hydraulic  gradient.  All  three  other  dimensionless  parameters  are  constant  for  simulations  2  and 
16.  It  can  be  seen  that  for  the  case  of  a  higher  hydraulic  gradient  (Run  16),  the  absolute  flow  rate 
through  the  gate  (Qgate)  increases.  However,  the  relative  flow  though  the  gate,  {Qrel),  and  the 
relative  size  of  the  capture  zone  remain  essentially  unaffected  by  the  increase  in  hydraulic 
gradient. 

To  examine  the  effect  of  variations  in  the  — ratio,  several  comparisons  can  be  made. 

^aquifer 

Comparison  of  Runs  6  and  16,  7  and  14,  and  8  and  12  are  three  comparisons  where  all  other 
dimensionless  parameters  are  constant.  All  three  comparisons  show  that  for  higher  values  of 
Kgate  relative  to  K^quifer^  both  Qgate  Qrel  increase  as  well  as  the  relative  size  of  the 
capture  zone.  Higher  hydraulic  conductivities  tend  to  draw  more  flow  towards  the  gate,  thereby 
increasing  the  flow  rates  through  the  gate  (both  absolute  and  relative).  Moreover,  the  size  of  the 
capture  zone  is  increased. 

Comparison  of  Runs  2  and  12  and  4  and  16  allows  an  examination  of  the  effect  of  the  — 

ratio  (the  ratio  of  the  width  of  a  single  cutoff  wall  to  the  depth  of  the  funnel-and-gate  system). 
Higher  values  of  this  ratio  result  in  lower  Qy^i  and  Qgate  values  and  smaller  relative  capture 
zones.  This  is  because  larger  values  of  this  ratio  result  from  higher  wf  values,  and  lower  df 
values,  in  which  case,  each  cutoff  wall  is  wider  and  shallower,  thereby  increasing  the  component 
of  flow  that  is  diverted  under  the  cutoff  walls  rather  than  through  the  gate. 

Finally,  the  ratio - -  (i.e.,  the  ratio  of  the  total  cutoff  wall  width  to  the  width  of  the 

gate)  can  be  examined  by  comparing  Runs  5  and  12  and  3  and  16.  Higher  values  of  this  ratio 
result  in  higher  Qgate  values,  but  lower  Qj-^l  values  and  smaller  relative  capture  zones.  The 
reason  for  this  is  that  larger  values  of  this  ratio  result  in  more  cutoff  wall  material  per  unit  of  gate 
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opening.  This  will  lead  to  more  flow  directed  through  the  gate,  but  the  increase  in  flow  through 
the  gate  is  not  proportional  to  the  increase  in  cutoff  wall  material.  This  results  in  a  lower  Qrel 
value. 
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(a)  plan  view 


(b)  cross  section  of  a  partially  penetrating  system 
Figure  1 .  Schematic  Diagram  of  a  Typical  Funnel-and-Gate  System. 
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Figure  3.  Streamlines  in  the  Vicinity  of  Funnel-and-Gate  System.  The  View  is  Looking  in  the  Direction 

of  the  Flow,  Towards  the  Funnel-and-Gate.  The  Dots  Represent  the  Beginning  of  Each  Streamline. 
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Figure  7.  Streamlines  and  Capture  Zone  for  Case  3  of  Comparison  Simulations. 


gate 


X 

=  widtli  of  one  (of  two)  cutoff  walls 
Wg  =  width  of  the  gate 
df  =  depth  of  the  funncl-and-gate  system 

Figure  9.  Schematic  of  an  Example  of  the  Funnel-and-Gate  Parameters  (Looking  in  the 
Direction  of  Flow,  Towards  the  Funnel-and-Gate). 
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Figure  10.  Streamlines  and  Capture  Zone  for  Rxm  #1  of  Run  Matrix. 
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Figure  12.  Streamlines  and  Capture  Zone  for  Run  #3  of  Run  Matrix. 
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Matrix. 


Figure  15.  Streamlines  and  Capture  Zone  for  Run  #6  of  Rim  Matrix. 
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Figure  16.  Streamlines  and  Capture  Zone  for  Run  #7  of  Run  Matrix. 


Streamlines  and  Capture  Zone  for  Rxm  #8  of  Run  Matrix 


Figure  18.  Streamlines  and  Capture  Zone  for  Run  #9  of  Run  Matrix. 
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Figure  19.  Streamlines  and  Capture  Zone  for  Run  #10  of  Run  Matrix. 
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Figure  20.  Streamlines  and  Capture  Zone  for  Run  #1 1  of  Run  Matrix. 


Figure  21.  Streamlines  and  Capture  Zone  for  Run  #12  of  Run  Matrix, 
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Figure  22.  Streamlines  and  Capture  Zone  for  Rim  #13  of  Run  Matrix. 
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Figure  23.  Streamlines  and  Capture  Zone  for  Run  #14  of  Run  Matrix. 
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Figure  24.  Streamlines  and  Capture  Zone  for  Run  #15  of  Run  Matrix. 


Figure  25.  Streamlines  and  Capture  Zone  for  Run  #16  of  Run  Matrix, 


APPENDIX  A 


SUPPLEMENTARY  OUTPUT  -  CASE  1  (OF  COMPARISON  CASES) 
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Figure  A-1.  XY  Slice  of  Hydraulic  Head  and  Velocities,  Z  =  9.5  Meters. 
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Figure  A-2.  XY  Slice  of  Hydraulic  Head  and  Velocities,  Z  =  5.0  Meters. 
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Figure  A-3.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  50.0  Meters  (Through  Gate). 


_J 

o 

4-> 


DZ 


LH 

tS 

(S 

G 

S 

>> 

J3 

in 

CN 

a 

Q 


o 


*a 

o  Q 

QJ  o 

31  S 
--  S 
u 

■  —  in  ^ 

Zi  o 

o  o 
+-’  I 

-a  c  a 

O  s- 

X  CJ  o 


C  -M  -M 
O  D  QJ 
•  -  un 

+-'  QJ 

o  u  o 

c  4-^  -  -  4-’ 
3  QJ  —  D 

LU  LH  (_n  Q 


.® 

LO 


tin 

CM 


44 


Figure  A-4.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  44.0  Meters  (Through  Cutoff  Wall). 
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Figure  A-5.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  40.0  Meters  (Through  CutotF  Wall). 
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Figure  B-3.  XY  Slice  of  Hydraulic  Head  and  Velocities,  Z  =  9.5  Meters. 
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Figure  B-4.  XY  Slice  of  Hydraulic  Head  and  Velocities,  Z  =  5.0  Meters. 
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Figure  B-5.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  50.0  Meters  (Through  Gate). 
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Figure  B-6.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  44.0  Meters  (Through  Cutoff  Wall). 
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Figure  B-7.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  40.0  Meters  (Through  Cutoff  Wall). 
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Figure  C-1.  XY  Slice  of  Hydraulic  Head  and  Velocities,  Z  =  9.5  Meters. 
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Figure  C-2.  XY  Slice  of  Hydraulic  Head  and  Velocities,  Z  =  5.0  Meters. 
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Figure  C-3.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  50.0  Meters  (Through  Gate). 
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Figure  C-4.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  47.0  Meters  (Through  Cutoff  Wall). 
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Figure  C-5.  YZ  Slice  of  Hydraulic  Head  and  Velocities,  X  =  44.0  Meters  (Through  Cutoff  Wall). 


